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Abstract 

Starting from an electrical dissipative illuminated one-diode solar cell with a given model data at room temperature 
(IsoVoo l^so: Fisho: Imax)': we present under physical considerations a specific mathematical method (using the Lambert 
function) for unique determination of the intrinsic electrical parameters (n, /p/,, R5, R,/,). This work proves that for a 
given arbitrary fixed shunt resistance R5/,, the saturation current Is and the ideality factor n are uniquely determined 
as a function of the photocurrent /p/,, and the series resistance R^. The correspondence under the cited physical 
considerations: R5 does not exceed ]0, 20[O and n is between ]0, 3[ and Iph and 4 are arbitrary positive {lR\), is biunivocal. 
This study concludes that for both considered solar cells, the five intrinsic electrical parameters that were determined 
numerically are unique. 

Keywords: Solar cell model; Electrical parameters; Electrical characterization; Lambert function; Shokley's equation; 
Numerical modeling 



Introduction 

Although the electrical dissipative one diode model has 
a potential of improvement in the efficiency and the sta- 
bility of the solar cell structure under illumination, to 
our knowledge the uniqueness and the authenticity of 
the extracted intrinsic electrical parameters associated to 
the model have not been studied previously. 

In this work we attempt to develop this concept and prove 
the uniqueness of the determination of these parameters. 

The one-diode model gives sufficient efficiency for 
earthly applications (Charles 1984). A precise numerical 
method using this model was presented in the early 
1980s by Charles et al. (1981; 1985). 

The use of the Lambert W-Function proposed by Corless 
et al. (1996) allowed demonstrating explicitly the Shokley's 
modified eq. (1) which is related to the equivalent elec- 
trical circuit model as shown in Figure 1. 



Where 7^/, is the photocurrent, n is the diode ideality 
factor of the junction, /jis the reverse saturation current, 
7?sis the series resistance and R^h is the shunt resistance. 

Each of these parameters is connected to the suited in- 
ternal physical mechanism acting within the solar cell. 
Their knowledge is therefore important. 

Several methods were proposed to determine the in- 
trinsic electrical parameters: Ipp,; n; 1^, R^; Rsh presented 
in eq. (1) of the solar cell. In particular, Jain and Kapoor 

(2005) established a practical method to determine the 
diode ideality factor of the solar cell. 

Ortiz-Conde et al. (2006) have used a co-content func- 
tion to determine these parameters. Jain et al. (2006) de- 
termine these parameters on solar panels. Chegaar et al. 

(2006) have used four comparative methods to deter- 
mine these parameters. 

More recently, Kim and Choi (2010) have used an- 
other method to determine the intrinsic parameters of 
the cell by making a remarkable initialization of the 
ideality factor n and the saturation current Is (Kim & 
Choi 2010). 
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Theoretical study: problem formulation 

To determine the solar cell intrinsic electrical parame- 
ters (k. Is, Iph, Rs, Rsh)> we put together a system of five 
equations (Lemma 2), and, solved by two different 
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numerical methods. The Lambert W-function is the re- 
verse of the function F defined from C+ in C hy F (W) = 
for every W in C+. 
Lemma 1: The Lambert W-function was derived from 
eq. (1) by expressing the current / in function of the 
voltage Vand vice-versa, as follows 



W 



V 



nVT(R,+Rsh) 



Rs + Rsh 



Rsh{ls 



Rs + R: 



■sh 



V = -I{Rs + Rsh) + Rshlph 
^IsRsh exp 



fro f Rs,{-I+I.+I,,) \\ 



-W 



V 



nVx 



nVj 



nVj 



(2) 



hVt +Rshh 



(3) 



p{i)=i 



\ 



\ 



nVj- 



\ 



nVr+Rshls 



J 



(4) 

We consider the following I (V) solar cell characteris- 
tics under illumination in generator convention as pre- 
sented in Figure 2. 

Where I^c and represent the short-circuit current 
and the open-circuit voltage respectively, R^ho is the 
slope of the I-V curve at the (0, Isc) point, Rso is the slope 
of the I-V curve at the (Voc, 0) point and I^ax is the max- 
imum power current, and Iph, 4, n, R^, and Rsh are the in- 
trinsic electrical parameters that should be determined. 

In order to simplify the problem formulation, we adopt 
the following abbreviations 

^~{Isc, Vac i-'^sO i^sdO /max ) ! {n,Is J ph fis '^sh ) 



Aj—exp 
A3— exp 



Rsh Rs (Is+Iph) 
n Vt {Rs+Rsh) 

Rsh {-Isc+h+Iph 
n Vt 



A4—exp 



Rsh (Is+Iph) 

n Vt 

Rsh (- ma 



^Is+I, 



pit, 



n Vt 
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From eq. (2) and at the point (0, Isc) we obtained 



Rs 



(5) 

Idem from eq. (3) and at the point (Vqc, 0) we obtained 

f2{X,Y)=Rsh Iph-W 'I n Vt^R^h h-Voc 
\ n v J J 

(6) 

The slope at the point (Vqc, 0) of the eq. (2) we obtained 



f^{X,Y)=-R,h- 



-Rso+Rs+Rsh 



1+W (^) 
The slope at the point (0, Isc) of the eq. (2) gives 



(7) 



f4{X,Y)=-R 



■sh ' 



Rsho+Rs+R 



■sh 



(8) 



1+W (^^) 
For differentiating eq. (4) and at the point (I = Imax) stems 



I 



W 



-Rs-Rsh+Rsh- 



Is Rsh \ 

n Vt 



\ 



1+W 

~^^max {Rs~^Rsh ) ~^Rsb^ph 



h Rsh ^4 

n Vt 



(9) 



Lemma 2: We have the following system 

f/,(x,r)=o 

f3{X.Y)=0 
fsi^,Y)=0 



(10) 



Proof: For / = and V=0, eq. (2) implies that/j(X, Y) = 
0 and for V= and / = 0 eq. (3) implies that/sCX, Y) = 0. 

The differential resistances: Rso and Rsho lead to the 
following two equations: fsiX, Y) = 0 and f4{X, Y) = 0. 

From eq. (4), maximal power obtained by: (Ij);-!,^^ = 0 
implies that/sK Y) = 0. 

In order to solve the system presented in Lemma 2 
(eq. 10) and determine the intrinsic electrical parameters, 
a set of experimental measurements (data) were used 
(Table 1). 

These measurements were collected from two dif- 
ferent solar cells under AMI illumination (E = 1 S = 
100 mW/cm^) at room temperature. 



Table 1 SAT and Cu2S-CdS cells experimental data 

Experimental data SAT cell (E = 1 S) CusS-CdS cell (E = 1 S) 



Vac (V) 


0.536 


0469 




0.45 


6.857 


L (A) 


0.1025 


0.04075 


^max (A) 


0.0925 


0.025 


RshO (^) 


1000 


41.905 


Vt(V) 


0.025875 


0.023527 



Our study concerns p-n junctions at both homo- and 
hetero-junctions: For the homo-junction, a 4 cm^ blue 
type monocrystalline silicon cell produced by SAT 
(1980) was used. For the hetero-junction we have used a 
frontwall Cu2S-CdS cell produced by a wet (Cleveite) 
process with significant losses of 4.28 cm^ square area. 
Two different numerical methods were applied in order 
to prove their authenticity. 

Numerical approach of the intrinsic parameters 

Newton's method 

The following function was considered 

F(X, Y) = (/i(X, Y),MX, Y)J^{X, Y),UX, Y)J,{X, Y)). 
Let Jf denote the Jacobian matrix defined by 



l<i,/<5 



So, Newton's method can be formulated as follows: For 



,0 jO jO 



p[,,Rs,Rsh^ as an initial condition and for 

all = 0, 1 ... until convergence; we have to resolve the un- 
known variable }* using the following system of equa- 
tions: Jp (Y'^) 5 Y'^ = - F (Y'^), where: + ^ = Y'^ + 5 Y'^ and: 

(nM^I^„Rk,Rk, 

In order to apply the Newton's method to this system an 
iterative program was developed in a MAPLE environment 
Monagan et al. (2003) using an accuracy of 20-digits. 

It depends on the choice of the initial data Y" by mak- 
ing sure that Jf (i*) 0 and by continuing the iteration 
process until a quadratic convergence is reached. 

At each increment, the program performs a test be- 
tween two successive iterations by assessing the Euclid- 
ean norm of their difference. The program was designed 
to stop the calculation when the test reaches a value 
smaller than the pre-set tolerance value. 

Hooke-Jeeves's method 

The Hooke-Jeeves method is based on numerical calcu- 
lation of the minimum of a function G without the use 
of gradient. This method is widely used in applications 
with convex G. 
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This method was used in this study to find the zero of 
function G (eq. 11) by minimizing in X and ysuch that 

GiX,Y)=0; G {X,Y)= \fAX,Y)\+ \f2{X,Y)\ 

MfsiXJ) I + \f4iX,Y)\+ \fs{X,Y) I 

(11) 

We recall that G (X, Y) = 0 is equivalent to the system 
presented in eq. (9) which leads to the determination of 
the intrinsic electrical parameters. 

This method has the advantage of being easily pro- 
grammed except the need to calculate gradient G. 

Existence and uniqueness of the solution 

To determine the existence and the uniqueness of the sys- 



implicit functions theorem where: H represents a continu- 
ously differentiable real-valued functions defined on a do- 
mam D in IR^x IR^ into 



hiylph, Rs-n, ls)=-Isc+lph-— 



Rsh 



exp 



nkT 



and 



hilph, Rs, «, h)=lph--jp^-h 



exp 



nkT 



By using the following notations: A = [Iph, Rs); B = («, Q 
Let /^(A, B) be the following Jacobian matrix: 



l<i,;<2 



tern presented in 


lemma 2 (eq. 9), 


we use the following 


Table 3 Frontwall CujS-CdS solar cell's intrinsic electrical 


Table 2 SAT solar cell's intrinsic electrical parameters (E = 1 S) 


parameters (E = 1 S) 






Intrinsic parameters 


Newton's method 


Hooke-Jeeve's method 


Intrinsic parameters 


Newton's method 


Hooke-Jeeve's method 


Iph (A) 


0.102502 


0.102002 


Iph (A) 


0.045528 


0.045487 


4 (A) 


5.987171985 X 10"' 


5.97501 X 10"' 


Is (A) 


8.2455 X 1 0"*^ 


8.0 X 10"*^ 


n 


1 .709464 


1.721481 


n 


2.183476 


2.176611 




0.016437 


0.016437 




5.355408 


5.298477188 


Rsh (O) 


1014.244754 


1000.412260 


Rsh (O) 


49.838828 


49.175974 
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Table 4 SAT solar cell's calculated I 
method 



-V values by Hooke's 



1/(1/; 

0.000000 
0.1 00000 
0.150000 
0.200000 
0.250000 
0.300000 
0.325000 
0.350000 
0.375000 
0.400000 
0.425000 
0.450000 
0475000 
0.500000 
0.536000 



/exp (A) 

0.102500 
0.102500 
0.102500 
0.102500 
0.102500 
0.101500 
0.101200 
0.100500 
0.099500 
0.09770 
0.09450 
0.08900 
0.07780 
0.05750 
0.00000 



'wookc (A) 

0.102501 
0.102396 
0.102333 
0.102245 
0.102079 
0.101669 
0.101241 
0.100510 
0.099246 
0.097048 
0.093216 
0.086533 
0.074895 
0.054718 
0.000000 



D (%) 

0.000009 
0.001015 
0.001631 
0.002494 
0.004124 
0.001665 
0.000405 
0.000099 
0.002559 
0.006718 
0.013774 
0.028509 
0.038787 
0.050842 
0.000000 



Let: {A°, B°) be a point in D such that H {A°, B°) = 0, 
and jfi{A°, B°) is invertible i.e. {Det (j^ (A°, B°));^0). 

The last step is to determine the neighborhood U x V 
where the following determinant of the Jacobian matrix 
will remain 



Det {jf,iA,B)) = 



-RJsc exp 



Voc exp 



Vqc 



exp 



nVj 



nVj 



1- exp 



nVr) 



Table 5 SAT solar cell's calculated l-V values by Newton's 
method 



V(V) 



le.n (A) 



'Newton 



(A) 



D (%) 



0.000000 
0.100000 
0.150000 
0.200000 
0.250000 
0.300000 
0.325000 
0.350000 
0.375000 
0.400000 
0.425000 
0.450000 
0475000 
0.500000 
0.536000 



0.102500 
0.102500 
0.102500 
0.102500 
0.102500 
0.101500 
0.101200 
0.100500 
0.099500 
0.09770 
0.09450 
0.08900 
0.07780 
0.05750 
0.00000 



0.102500 
0.102397 
0.102337 
0.102255 
0.102104 
0.101739 
0.101359 
0.100708 
0.099580 
0.097614 
0.0941 73 
0.088149 
0.077613 
0.059257 
0.000000 



0.000000 
0.001005 
0.001592 
0.002395 
0.003878 
0.002354 
0.001571 
0.002069 
0.000804 



0.003472 
0.009654 
0.002409 
0.030556 
0.000000 



Table 6 Frontwall CujS-CdS solar cell's calculated 1- 
by Hooke's method 


-V values 


V(V) 


W (A) 


'wooke (A) 


D(%) 


0.000000 


0.0408000 


0.041010 


0.0051470 


0.050000 


0.0393 


0.039746 


0.01 1348 


0.100000 


0.0373 


0.038120 


0.021983 


0.200000 


0.0315 


0.032689 


0.037746 


0.250000 


0.0273 


0.028445 


0.041 941 


0.275000 


0.0250 


0.025929 


0.037160 


0.300000 


0.0225 


0.023173 


0.02991 1 


0.325000 


0.0196 


0.020201 


0.030663 


0.350000 


0.0165 


0.017038 


0.032606 


0.375000 


0.0132 


0.013705 


0.038257 


0400000 


0.0099 


0.010225 


0.032828 


0450000 


0.0025 


0.002895 


0.158000 


0469000 


0.000000 


0.000000 


0.000000 


This determinant does not depend on Ipf, and is linear 
with /j. The Rs and n dependences of the determinant 
are illustrated in the following figure (Figure 3). 

The minimum of the determinant in ] 0, 20[x]0, 3 [ is 
10'^. Consequently the investigated neighborhood U x 
Vis IR;x]0,20[x]0,3[xIR;. 

The implicit functions theorem gives the existence of a 
unique function B = (p (A) defined in U into V of class 
and for any {A,B)gUxV, H {A, (p {A)) = 0. As a 
result the (p Jacobian matrix is given by the formula: 
(A) = Jfj (A, (p (A))"V^ (A, (p (A)) and consequendy, we 
prove for a given arbitrary fixed shunt resistance R^h, that 
the saturation current 1^ and the ideality factor n are 
uniquely determined in function of the photocurrent Ipf,, 
and the series resistance R^- 


Table 7 Frontwall Cu2S-CdS solar cell's calculated 
values by Newton's method 


l-V 


V(V) 


Up (A) 


^Newton (^} 


D(%) 


0.000000 


0.0408000 


0.040750 


0.001226 


0.050000 


0.0393 


0.039420 


0.003053 


0.100000 


0.0373 


0.037667 


0.009839 


0.200000 


0.0315 


0.031870 


0.01 1 746 


0.250000 


0.0273 


0.027522 


0.008131 


0.275000 


0.0250 


0.025000 


0.000000 


0.300000 


0.0225 


0.022273 


0.010191 


0.325000 


0.0196 


0.019362 


0.012292 


0.350000 


0.0165 


0.016290 


0.012891 


0.375000 


0.0132 


0.013075 


0.009560 


0400000 


0.0099 


0.009737 


0.01 6740 


0450000 


0.0025 


0.002748 


0.099200 


0469000 


0.000000 


0.000000 


0.000000 
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Experimental and theoretical results, discussion of To prove the authenticity of the model, we should 

related authenticity calculate the current / listed as Ith by the use of the 

Tables 2 and 3 list the intrinsic electrical parameters obtained intrinsic parameters at different points 

values of the two cells determined by Newton's of the I-V curves. These points are compared 

method and Hooke-Jeeves's. with the corresponding experimental current values 



0 12 
008 

< 

006 
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0<M 



0 02 
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0.045 
0 

0.035 
0.03 
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r 0.025 h 
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0.015 
0.01 

0.005 
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Voltage : V (V) 




06 
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- • — theoretical characteristic by Xe\«on method 




Cii2S-Cd.S Solar CeU 
E=1S 



0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45' 0.5 

Voltage :V (V) 

Figure 4 Experimental I-V Characteristics. (A) c-Si blue SAT solar cell. (B) Frontwall Cu2S-CdS solar cell. 
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listed as lexp- The accuracy is evaluated by the par- between experimental and theoretical results are also 
ameter D (%). The values of the called accuracy D listed in Tables 4, 5, 6 and 7 and does not exceed 
{%) corresponding to the percentage deviation 0.2%. 



X 10" 



< 




Hooke exp 
'Newton 'exp' 



0.3 0.4 
Voltage : V ( V) 



0.7 



B 



X 10" 



Hooke exp 
Newton exp 



Cu2S-CdS Solar Cell E=1S 
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1 

0.5 

O' 

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Voltage : V ( V) 

Figure 5 Absolute error between experimental and calculated current. (A) c-Si blue SAT solar cell. (B) Frontwall Cu2S CdS solar eel 



o 
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The graphs presented in Figure 4A and B show how 
close the values calculated by the two used numerical 
methods to the experimental ones. 

These figures are very sensitive to the effects of the circuit 
parameters with localized constants and especially to the 
quality of the cell. Figure 5A and B outline the absolute er- 
rors between the experimental and calculated current as a 
function of the cell bias voltage by the two numerical 
methods. Although D values of the SAT solar cell in the 
state of the art are weaker than those of Cu2S-CdS solar cell 
with significant losses; absolute error (Figure 5A) goes to a 
maximum at Voc-neighborhood. This maximum is weaker 
in the case of Newton's method, so denoting a better con- 
vergence of this method compared to Hooke's. Although in 
the case of the Cu2S-CdS solar cell with significant losses 
this indeterminacy on Rs disappears, the calculated TV 
curves show a better convergence of Newton's method. 

Conclusion 

In this study a simple and specific method (without ap- 
proximations) was proposed to extract intrinsic electrical 
parameters of the one-diode solar cell model under 
AMI illumination (IS). 

The proposed approach includes parasite and dissipa- 
tive elements such as series resistance Rs and shunt re- 
sistance Rsk- 

The use of the Lambert W-function has allowed to ex- 
press explicitly the current / as a function of the voltage 
V^from the modified Shockley's eq. (1). 

However, it is important to highlight that the proposed 
method is valid for all measured I-V characteristics 
under any illumination intensity. 

The implicit functions theorem was used to demon- 
strate the uniqueness of the solution. The physical consid- 
erations of the problem have also been taken into account. 
This procedure has proved the uniqueness of the solution. 

Two different numerical methods: Newton's method 
and Hooke-Jeeves's were used to determine these param- 
eters and reconfirm the uniqueness of the solution. 

To prove the authenticity of this extraction method, two 
different types of solar cell structure were used: a SAT mono- 
crystalline silicon homostructure in the state of the art, and a 
frontwaU Cu2S-CdS heterostructure with significant losses. 

Moreover, as MATLAB has limitations toward large 
numbers manipulation (> exp (100)), MAPLE software 
was selected for this calculation. 

For the two cell types, both used numerical methods 
converge in each of cases, towards two series of theoretical 
results with relative accuracy about 3% in the case of the 
weak series resistance. 

Nomenclature 

T: Thermodynamic Temperature in Kelvin (K) 
q: Electron Charge = 1.602*10'^'' C 



k: Boltzmann constant = 1.38*10'^^ J/K 

Vji Thermal voltage = kT/ q 

Voc- Open circuit voltage 

Isc'- Short-current voltage 

/: Output current 

V: Output voltage 

imax'- Maximum power current 

Vmax- Maximum power voltage 

Pmax- maximum power 

Ipk- Photocurrent 

Is- Diode reverse saturation current 

loc- Calculated current at the (Voc, 0) point 

Vsc- Calculated voltage at the (0, Isc) point 

n: Diode quality factor 

Rsh- shunt resistance 

Rsho- Differential Resistance at the (0, I^c) point 
Rs- Series resistance 

Rso- Differential resistance at the (Voc, 0) point 
W: Lambert's function 

C+: the set of complex numbers with positive real part. 
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